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Abstract 

A composition of an integer n is called Carlitz if adjacent parts are dif- 
ferent. Several characteristics of random Carlitz compositions have been 
studied recently by Knopfmacher and Prodinger. We will complement 
their work by establishing asymptotics of the average number of distinct 
part sizes in a random Carlitz composition. 







1 Introduction 



In this note we obtain precise asymptotics, as n ^ oo, for the expected number 
of distinct part sizes in a random Carhtz composition of an integer n. Let us 
recall that a tuple (71, . . . ,7^) is a composition of an integer n if the 7j's are 
positive integers, called parts, such that 7^ = n. The number k is the num- 
ber of parts and the values of 7j's are called part sizes. There are 2"'"^ different 
compositions of n. A composition is called Carlitz if the adjacent parts are 
different, i.e. if 7j 7^ Ij+i for j = l,...,fc — 1. For example, out of sixteen 
compositions of the integer 5, seven are Carlitz, namely (5), (4, 1), (1,4), (3,2), 
(2,3), (1,3,1), and (2,1,2). For other values, see ||, Sequence A003242]. We 
denote the set of all Carlitz compositions of n by f2„. Carlitz compositions have 
been introduced by Carlitz ||l| who found the generating function for the total 
number of them. They have been subsequently studied by Knopfmacher and 
Prodinger ||] (see also |^). These authors found the asymptotics of the total 
number of Carlitz compositions. They also studied several parameters (like the 
number of parts, the size of the largest part among other things) for random 
Carlitz compositions. "Random Carlitz composition" means a composition cho- 
sen according to the uniform probability measure on f2„. This measure will be 
denoted by P and E will denote integration with respect to P. In this set- 
ting, various parameters of Carlitz compositions become random variables and 
their probabilistic properties are to be studied. For example Knopfmacher and 
Prodinger found, among other things, the exact asymptotic behavior of the ex- 
pected number of parts and the expected size of the largest part. One question 
that they left open concerned the expected value of the number of distinct part 
sizes. Dm and the purpose of this note is to answer their question. In order 
to state the result we need to introduce some more notation and we will try to 
closely follow the notation of Knopfmacher and Prodinger. First, the number 
of distinct part sizes is defined formally as follows: if (71, . . . , 7^) is in fJ,, then 

k 

Dn = 1 + 3 = 1,. ..A-1}, 

1=2 

where Ia denotes, as usual, the indicator function of the set A. Secondly, we 
introduce a function ct of a complex variable defined by 

00 j 

The equation a{z) — 1 has the unique real solution p = 0.571349 ... on the 
interval [0,1]. (The relevance of this is that, as was shown by Carlitz, the 
generating function of Carlitz compositions is equal to 1/(1 — a{z)).) We have 

Theorem 1 With the above notation, and letting { ■ } denote the fractional 
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part, as n —> oo, we have 

\n{n/a'{p)) 1 7 , u i! ^ t / V M^ , n^ 
^ ln(l/p) + 2 - MTM + + ''^l)' 

where 7 is Euler's constant, ho is a mean zero function of period 1 whose Fourier 
coefficients are given by 

1 , 27rie , 

Approximating all the constants and using the fact that the gamma function 
decays very fast along the imaginary axis, it can be seen in particular, that 

ED„ = Ci Inn - C72 + hoi{pHn/a' (p)}) + o(l), 

where Ci — 1.786495 . . ., C2 — 2.932545 . . ., and the amplitude of Hq is bounded 
by 0.5882304 •• • x 10"^. 

The approach is as in |^ via generating functions. We let Ij denote the set 
of those Carlitz compositions that contain at least one part of size "j" . Using, 
without any risk of confusion, the same notation for a set and its indicator we 
have 



and therefore 



where A"^ denotes the complement of a set A. Denoting by a„ and a„ j- the num- 
ber of all Carlitz compositions of n and the number of those Carlitz compositions 
of n that do not use size j we have 

P(/n = 

an 

The sequence (a„) was studied in Knopfmacher and Prodinger so we need 
only to study the numbers a„j . In order to do that we will build their generating 
function. The construction follows the ideas of but since that paper is short 
on some details and may be a bit difficult to read for a new adept, we will 
provide a fairly detailed argument. 

We would like to mention that the method of our paper does not appear 
to be strong enough to yield information about the limiting distribution of 
the variable A bivariate generating function for Z)„ would be extremely 
welcome. But, we were unable to get it. Although it is known (c.f. |^) that the 
total number of parts in Carlitz composition satisfies the central limit theorem, 
we think that it is unlikely that D„ will have the same property. It would be 
very interesting and desirable to find the limiting distribution of Z?„. 
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2 Generating function 

In this section we will prove the following statement 

Proposition 2 Let Cj {z) be the generating function of the sequence {unj , n > 
0}. Then, 

1 - CTj(z) 

where 



e>i 



Proof: Let a^*^^ (n, m) denote the number of Carlitz compositions of n with the 

properties: 

• they have exactly k parts 

• they do not contain a part of size j 

• the last part is of size m. 
Then for > 1 and m ^ j we have 

^(/s+i) _ ^(fe) _ ^-j _ — m,m), 

where a^'^'' (£) is the number of Carlitz compositions of i into k parts, none of 
them equal to j. Since j is fixed throughout the argument, for the ease of 
notation we supress the subscript j throughout the argument. We additionally 
require that aP'\t) and a^''^i,m) vanish whenever i <0. Let 



fk{z,u)= a^'"\n,m)z'' 



m>l, n>0 



Since the compositions enumerated in a^^\n,m) do not contain a part of size 
j, we have a^''\n,j) = 0. Hence 



fk+i{z,u) = 5^ 5^a('=+i)(n,m)^V 

ra>Om>l 

„>0 m>l „>0 

- m^j - 

= E E (a('=)(n-TO)-a('=Hn-m,m))2"u'" 

n>0 ">>i 

= E E a('=)(n-m)z"u'" 

n>0 ">>! 

- E E a^''\n - m, m)z''u'^ 



n>0 m>i 
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The first sum above is equal to 

n>0, m>l n>0 

t>Q, m>l m>l 
i — 2:m 

where, in the second step we changed the summation indices n — m = £ and 
n — j = m, respectively. By the same argument 

Y a^''\i,m)z'{zur = Mz,zu). 
m>l, e>o 

Thus, for A; > 1 we have 

(zu ■ \ 
{zuy - fk{z,zu). 
1 — zu J 

Letting /o(-2, u) = \ the last line can be rewritten as 

(ZXl ■ \ 
[zuY - fk{z, zu) + 4,0, 
1 — zu J 

for > 0. Introducing the function 

fe>i 

and summing over > we obtain 

F{z, u) = F{z, 1) f - {zuy) + (-^ - {zuy) - F{z, zu). 



^1 — zu ) \\ — zu 

This equation can be iterated to yield 



Hence, for u = 1 we get 

F{z,l) = F{z,l)<Tj{z) + aj{z), 



where (Jj{z) is defined by (^. Finally, letting 

Q(z) = l + F(z,l), 

we see that 

C,(z)-l = (Q(z)-l)a,(z) + a,(z), 

which, since 

F{z, 1) = ^ 1) = ^ ^ a^'''\n,m)z^^ = ^ a„,,z", 

fc>l fc>l ">o n>0 

~ ~ m>l ~ 

means that 

CA^) = 1 \ V 

l-crj(z) 

is the generating function of the sequence (anj)- 

3 Singularities of the generating function 

A starting point of our analysis is the fact that the generating function of Carlitz 
compositions has the unique singularity in the disc {z : \z\ < 0.663}. This 
singularity is the unique real root, p, of the equation <7{z) = 1 on [0,1]. The 
numerical approximation of that root is p = 0.571349 . . .. Since a{z) and <Jj(z) = 
(t(z) — z^ /{I — z^), do not differ by too much, functions aj will have the same 
feature, at least for j's sufficiently large. In fact, on the disc {z : \z\ < 0.663} 
the functions satisfy the following: there exists 6 > such that for every j > 6 
the equation aj(z) = 1 has the unique real simple root on [0, 1]. Furthermore, 
these roots, which will be denoted by pj, have the following properties 

1. Vj >6,0<pj <p + S, 

2. Vj > 6 all roots ^ of (Jj{z) — 1 other than pj satisfy |^| > p + 26. 

3. Pj are strictly decreasing for j > 6 and pj p as j ^ oo. 

A justification as well as a discussion of the remaining case 1 < j < 5 is post- 
poned until the appendix. We will need asymptotics of pj and wc will use the 
"bootstrapping method" . Rewriting 

aj{z) = 1 

as 

a{z) = (3) 

letting Pj — p + Sj, where Sj = o(l) and substituting the latter expression for 
Pj into Q we get 

1 1 

a[p + £j 
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Using Taylor's expansion on both sides we see that 

Since a(^p) = 1 we obtain that Sj ~ p' /a'{p) + o{p'). Hence 

Pj- + (4) 

which is sufficient for our purpose. Let Aj = —l/<Tj{pj) be the residue of 
1/(1 - aj{z)) at pj. Then 

1 



l-crj(z) Z-pj 

is analytic in the disc {z : \z\ < p + 6} and by the Cauchy integral formula we 
get 

A,, „/ 1 



Finally, since 



aj{z) ^ (7{z) 



1-zJ' 
we get 

a',{p,)=a'ip,)^-^^^. 

(1 - p^jY 

Hence, taking into account we get 
Consequently, for j > 6 

for some (5 > (universal for j > 6) . 

4 Asymptotics 

The following claim will account for most of the asymptotic analysis of (^ 
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Lemma 3 A 



s n ^ (X, 



■ 1 \ / ■ n 



This statement immediately implies our theorem since the asymptotic be- 
havior of the series on the right is known (see e.g. [Q, 0, Section 7.8 and 
references therein]). 

Proof of Lemma 3: The sequence (a„) was studied in ||] and one has 



1 

Combining this with (|^) , for j > 6 

^ i + opy) 

an {^l + p^-^^-p^y 
where A ^ S / p > 0. Consequently, 



- + 0((1 + A)-"), 



■E 



'i- , ,.,. +o((i+ A)-)' 



We will show that 

where we have set aj — lla'{p) + o(l). To this end, we will split this sum as 



where j„ will be chosen momentarily. Since aj \/a'{p) > we have 
as long as jn — > oo at any rate. The other term of the sum is bounded by 



^E^l+Z^^'^j)"" < Cjnexp{-n\n{l + p^"-^aj)). 
3=e 
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Since ln(l + x) > x — a;^/2, for x > we get the bound 

Ci„exp ^-nap^""^ ^1 - — < Cj„ exp(-cnp^""^). 
Choosing j„ ~ /x logj^ n we see that this expression is bounded by 

0(logn-e-^"'"'*) =o(l), 

whenever ^ <1. 

It remains to consider the sum 



j=5 

We first replace this sum by a more convenient one 

n-l 

5^(l-(l-p^a,r). 

i=5 

The difference is at most 



For a jo which will be chosen momentarily, we consider 

+ p^a 



i=5 

Since each summand is nonnegative, this sum can be upperbounded by 

jo / J \ n jo / „■ - 

< {jo + 1) exp (- ) < (io + 1) exp (-cp^^n) , 

for some absolute constant c. By choosing jo so that 

cp^°n > Klnn, 

i.e. 

ln((cn)/(Klnn)) 

- HVp) ' 
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we see that the sum up to jo is bounded by clogn/n'^. For the remaining range 
of j's we write 



~ 1 + p^a\ 1 + piaJ 

where we have used the inequaUty (a + 6)" — a" < nb{a + &)"~-^ vahd for 
nonnegative numbers a and b. Hence 

The choice of jo so that p^" = 8 (log n/n) is within the previous constraint, and 
for that choice we have 

Using the same argument we can show that a.^ can be replaced by its limit 
a = l/(j'{p), i.e. that we have 

n— 1 n— 1 

^ (1 - (1 - p>air) = ^ (1 - (1 - p>ar) + 0(1). 
Finally, the sum on the right can be increased to 

oo 

^(l-(l-p^ar), 

3=0 

since 

^(1-(1-pM") <an^p^"=o(l) 

and is thus neghgible. Also, for each fixed j, with 1 < j < 5, an,j/an = o(l) as 
we will indicate in the appendix. This means that 

5 



and proves the Lemma. 
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5 Appendix 



In order to show that Gj {z) = 1 has a unique real root pj and the existence of 
a ^ > satisfying the asserted properties we rewrite a{z) in a more convenient 
form 

m— 1 

which can be done by expanding z^ /{I — z^) into geometric series and inter- 
changing the order of summation. Now, dj {z) = 1 can be rewritten as 



- + 2;™ I- z^ 

m— 1 



-1 = 0. 



We want to use Rouche's theorem. To this end spht the left hand side as 
f{z)+9j{z) where 

6 m 
m=l 



It can be verified that 



min 1/(^)1 > 0.28. 

|2|=0.663 



Actually Maple suggests that a stronger claim is true, namely that minimum of 
\f{z)\ on that circle is attained at z = 0.663 and is 0.283467. We have not tried 
to prove it and thus we claim only the weaker statement, which can be verified 
by evaluating |/| at sufficiently many points and using its Lipshitz property: for 

\zi\ = \Z2\ =r, 



6 I jy-^ jYf^ I 

1/(^1)1-1/(^2)1 < \fizi)-fiz2)\<J2wi^'' 

m=l ^ 



6 m — 1 

mr 



'){l + z^ 

E inr " 
m=l ^ ' 

where C < (1 - r)-^ mr™"! < (1 - ^)"*- Since for \z\ = 0.663 

m=7 ' ' ' ' 

- (1 - 0.6637)(1 - 0.663) ^ 1 - 0.6636 - ' ' 



we see that 



|z|=0.663 |z|=0.663 
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uniformly for all j > 6- Thus / + gj and / have the same number of roots in 
the disc \z\ < 0.663. But / has a unique (and thus necessarily real) root on that 
disc (this again may be verified by Maple after converting it to a polynomial 
equation). The existence of S follows since the same argument can be repeated 
for a disc with a radius just a bit smaller. As for the monotonicity, it suffices 
to rewrite (Jj{z) = 1 as cr(z) — 1/(1 — z^) and notice that <y{z) is increasing as 
a function of real variable and that 

1 1 

< 



1 - Z-Jl 1 - 

whenever ji > j2 and < z < 1 is real. 

The same argument (with changed parameters) can be used to verify that 
for j — 2, 3, 4 and 5, on the disc \z\ < r, < r < 1, aj{z) = 1 has a unique real 
root pj > p. The hardest case is j = 2. We found that the splitting 

m—1 m— 26 

will do the job. On the circle \z\ = 0.8 one has 1/2(2) | > 0.06, \g{z)\ < 0.016. 
The polynomial resulting from multiplying /2 by the product of denominators 
is of degree 264 (after cancellations). Maple (somehow reluctantly) shows that 
the root closest to zero is real and is about 0.78397 (there is another real root 
of about 0.927122) and the next closest to zero roots are complex conjugate 
and have absolute value around 0.81914. For j ~ 3, 4, 5 one can get away with 
letting g{z) — J2m=ii ^'"/(l + 2™)j and choosing r — 0.75. 

The same method could presumably be used to force the argument for j = 
1. However, here matters would be computationally worse. Furthermore, the 
equation cri (z) = 1 does not have real solutions, and thus, its closest to zero root 
(if it exists) would have to come from a pair of complex conjugates. Luckily, for 
that case we can use a different argument based on the probabilistic approach 
used in |^ (see also 0). Let us briefly sketch it. Consider the set C„ of all 2"~^ 
compositions of an integer n and let Q be the uniform probability measure on 
C„. Since the restriction of such a measure to any subset is again the uniform 
measure on that subset we can view the uniform measure P on the set of all 
Carlitz compositions as a conditional measure obtained by restricting Q to Qn, 
i.e. 

P( • ) = Q( • \^n). 

Let Al be the set of all compositions of n that do not use part size 1 and recall 
that /f is the set of all Carlitz compositions with this property. Then, 

an -P(A)-Q(A|iM- ^Q(f7„)- 
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Now, by the result of Knopfmacher and Prodinger 

Qi^n) > c-^ = c(0.875)", 

and we need to upper bound Q{A1). To this end we wiU use the observation 
made in or [Q that a random composition of n is distributed hke 

(fi,...,f,), 

where, with (F^) being a sequence of i.i.d. geometric random variable with 
parameter 1/2, we let 

fc 

r = infjfc ; Fg > n}, 

9=1 

and 

Tq if q < T, 



J2lJi Tg if q 



Furthermore, r is distributed like 1 plus a binomial random variable with pa- 
rameters n — 1 and 1/2. Therefore, 

r r — 1 

Q(^D = Q(n{rfc>i})<Q(n{rfc>i}) 

k=l k=l 
T-1 

< Q(T<fco) + Q(n{rfc>i}n{T>fco}) 



fe=i 

ko 



< Q(r-Er<fco-ET) + Q(f|{Ffc>l}) 



fe=i 

< Q(|r-Er|>^-fco)+(l-i''" 

< 2exp ^-2-^^ +exp(-fcoln2) , 

where we have put fcg (n+l)/2 — and 0„ = a(n+l), for some < a < 1/2. 
Then we get 

QiAl) < C (e-'"')" + C (^e-(i/2-")"^2j" 

To be able to claim that Q(y4!j;)/Q(ri„) tends to zero, we need both of the 
inequalities 

e"2"' < 0.875 and e^'^/^^^^'^^ < 0.875 
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to be satisfied simultaneously. But this can be achieved by choosing a to be 
any number subject to 0.26 < a < 0.3. This completes the argument. 
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Will, who brought the paper ||^ to our attention by asking the second named 
author whether the methods used in |^ could be used to rederive some of the 
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